^-measure and Disoriented chiral condensates 



B. Mohanty 

Institute of Physics, 751-005 Bhubaneswar, India 



Abstract 

Fluctuations in the ratio of neutral to charged pions arising due to formation of 
disoriented chiral condensates (DCC) are discussed using the ^-measure. The prop- 
erties of the measure for various cases of DCC and non-DCC are discussed. The 
effect of detector efficiencies and other experimental factors are presented. Applica- 
tion of the ^-measure to simulated data, within the context of a simple DCC model 
are also discussed. 



In heavy-ion collisions, there is rapid expansion of collision debris along the 
longitudinal direction, leading to supercooling of the interior interaction re- 
gion. This may result in the formation of domains of unconventionally oriented 
vacuum configurations as allowed by chiral symmetry called as disoriented chi- 
ral condensates (DCC) [1]. Detection of these would provide important infor- 
mations on the vacuum structure of the strong interaction and the nature of 
chiral phase transition. To have a system in which chiral symmetry is restored 
in the laboratory and to be able to study the above features is one of the main 
goals of ultra-relativistic heavy-ion collision experiments. 

It has been predicted that DCC formation is associated with large event-by- 
event fluctuations in the neutral to charged pion ratio. The probability of the 
neutral pion fraction [2,3], /, is 



P(f) = l/2Jf where f = N n o/N w , (1) 



N n o and N n being the total number of neutral pions and the total number 
of pions respectively. The corresponding distribution for non-DCC events is 
a gaussian with < / >= 1/3. It can easily be seen that for events with 
DCC there is a strong anti-correlation in the production of neutral to charged 
pions. Several theoretical calculations regarding various aspects of DCC exists 
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in literature, starting form the probability of DCC production [4] to life time 
of DCC [5]. Although there is an absence of a dynamical model of DCC, still 
the theoretical calculations are encouraging with respect to formation of DCC 
in heavy-ion collisions. 

A typical experiment, looking for DCC would consist of two detectors, one 
to detect charged pions and other to detect photons from the decay of 7r°'s. 
They must have a common r)-(f) overlap with rj coverage as much as possible. 
From the detected hit patterns one tries to see if there is any fluctuation in /, 
which would indicate presence of DCC-type fluctuations. Typical event struc- 
tures would be similar to the anti-Centauro events reported by the JACEE 
collaboration [6]. Results from other cosmic ray experiments have not ruled 
out the possibility of a DCC formation mechanism [7]. As far has accelerator 
based experiments are concerned, several experiments have attempted to look 
for DCC by colliding hadrons and heavy-ions. Hadron-hadron collision exper- 
iments like UA1 [8], UA5 [9], DO, CDF [10] and MINIMAX [11], having 
from 540 GeV to 1.8 TeV, have so far reported null results. Heavy- ion collision 
experiments like WA98 [12,13] and NA49 [14] at CERN SPS have so far put 
upper-limits on DCC production. In future several experiments have planned 
to look for this interesting phenomena at RHIC [15] and LHC [16]. 

Several techniques have been developed to look for DCC. These includes N 1 
vs. N ch correlation [13], "robust" variables [11] and those based on multi- 
resolution analysis techniques such as discrete wavelet transformations [13]. It 
must be mentioned that it is important to have the right variable or method 
to look for fluctuations which are exotic. The importance of this in general has 
already been emphasized by some of the recent calculations [17]. One of the 
important observables for looking at fluctuations is the ^-measure [18]. It has 
been developed specifically to remove the influence of trivial geometrical fluc- 
tuations and the effect of averaging over many particle sources. It's utility in 
looking for fluctuations in transverse momentum, azimuthal fluctuations [19] 
and studying chemical fluctuations [20,21] has been shown. More recently the 
usefulness of $ in studies of charge fluctuations have been demonstrated and 
it has been shown to have some advantages over other techniques [22] . Its lim- 
itations have also been discussed [23] . In this letter we try to see the usefulness 
of ^-measure for looking at DCC-type fluctuations. 

First we briefly recall some of the basic equations related to ^-measure in 
general and then study its properties for a DCC-type fluctuation. The $- 
measure for a system of particles is defined as, 



z is the single-particle variable, z = x — x, with x being the probability 
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(averaged over events and particles) that a produced particle is of the sort of 
interest, say it is 7r°. One easily observes that ~z — 0. While Z is the event 
variable, which is a multi-particle analog of z, defined as Z = f Y^=\{^i ~ x), 
where the summation runs over particles from a given event. By construction 
< Z >= 0, where < ... > represents averaging over events. 

As discussed in Ref. [20] one can compute $ for the system of particles of two 
sorts, a and b, e.g. n° and n ± . xi = 1 when i—th particle is of the a type and 
Xi = otherwise. The inclusive average of x and x 2 would then be given as 

x = Yl xP * = p i , ^ = E *? p x = p i . 

x=0,l i=0,l 

where Pi is the probability (averaged over particles and events) that a pro- 
duced particle is of the a sort. Thus, 

p t <iV " > 



<N a > + <N b > ' 



with iV a and N b being the numbers of particles a and 6, respectively, in a single 
event. One immediately finds that z = while 



5 - Pl " Pl - <iV>2 ' (3) 



where N = N a + N b is the multiplicity of all particles a and 6 in a single event. 
As the event variable Z is defined as N a — xN, one gets 

<Z> = <iV a >-Ic<iV>=0, 
< Z 2 > = < N 2 a > -2x < N a N >+x 2 <N 2 > . 

Which leads to 

< Z 2 >< N > 2 =< N b > 2 < N 2 > + < N a > 2 < N 2 > -2 <N a >< N b >< N a N b > , 
then 



< Z 2 > <N b > 2 . Ar2 2 \ ,,n 

TWT = Tn^ 2 {<n ° > - <n ° > ' (4) 

+^«iV l '>-< 1 V l >-) 

- 2 <jV °y f > < < ^ > - < jv. >< jv» > ) . 
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The fluctuation measure $ as given in Eqn. (2) is completely determined 
by Eqns. (3, 4). Now we identify N a as N n o and Nj, as N c h = N n + + N n -. 
If N n = N n o + N ch is the total pion, then we can write N n o = fN n and 
N ch = (1 — f)N w . Where / is the fraction of neutral pion out of total number 
of pions in a given event. With this we get, 



z 2 =< f > (1- < / >) (5) 
In order to calculate ^ 2 ^ , we consider the following equations. 

< SN 2 >=<f 2 >< 5N 2 >+<N 2 >< 5f 2 > (6) 
+2 < N n >< f >< 5N n 5f > 



< 5N 2 ch >=<{!- ff ><5Nl>+<Nl>< 5f > (7) 
-2 < N n >< (1 - /) >< 5N5f > 



< 5N w0 5N ch >=<f ><(!-/)>< 5Nl >-<N 2 x5f 2 >- (8) 
< N w >< f >< SNJf >+ < N w >< (1 - /) >< SNJf > 

where, < 5N >= N- < N > and < 5N 2 >=< N 2 > - < N > 2 . 

Using Eqn. 6, Eqn. 7 and Eqn. 8 in Eqn. 4 one can easily get the expression 

A small assumption is involved in the above derivation, that the fluctuations 
in N n is small and the term + ~ N n . 

Using Eqn. 5 and Eqn. 9 in Eqn. 2, the ^-measure can be defined as, 

$ = y/< ~N« >< 5f 2 > - ^< f > (1- < / >) (10) 

Let us now consider its properties within three simple models of the multi- 
plicity distribution. 
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Non-DCC case - It is known from pp experiments [24] that the produced 
pions have their charge states partitioned binomially with mean of / at 1/3. 
The fluctuation in / is inversely proportional to the total number of pions 
given by < 5f 2 > =< f > (1— < / >)/N w . Then one can easily see that, 

^non-DCC = (11) 

Non-DCC but N w o and N c h are correlated - If N w o and N C h are assumed to 
be correlated in such a way that there are no DCC-type fluctuations in the 
events. The neutral pion fraction, which is defined by the ratio A^o/A^, is 
assumed to be strictly independent of the event multiplicity. Then, f = a, 
where a is a constant smaller than unity. Then < 5f 2 > = and ^-measure 
is given as 

$ = -y/ a (l -a) . (12) 

DCC case - For the DCC case, the probability distribution is given by 
Eqn. (1). Using this, one can easily see that < / >= 1/3, as in the non-DCC 
case, while < 5f 2 >= 4/45. So the ^-measure is given as 

$ DCC = ^< N w > 4/45 - ^2/9 (13) 

For a typical case, where the total pion multiplicity, N n , in the experiment is 

300, we find for the DCC case the above ^-measure is about 4.7. 

This is a result where all the pions observed are of DCC origin. However a 

more realistic case [25] would be to include pions from non-DCC sources as 

well. 

Background pions - In a given event pions can originate both from a DCC 
source and from other non-DCC sources, even in a DCC event. In such a case 
the probability distribution for the neutral pion fraction, P(f) [26], is given 
as, 

P(f) — J df DC C dfnon-DCC P(fDCc) P(fnon-DCc) 

S(f - PfDCC " (1 " P)fnon-DCc) (14) 

where /3 corresponds to fraction of DCC pions out of a total number of pions 
N n , produced in the event, focc an d fnon-Dcc correspond to the neutral pion 
fraction for pions from DCC and non-DCC sources in the event, respectively. 
With the presence of non-DCC pions, one can see, the neutral pion fraction 
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Fig. 1. Variation of ^-measure as a function of the pion fraction (f3) from a DCC 
origin. The plot is obtained for an average pion multiplicity of 300. 



can no longer start from zero nor can reach the full value of unity, the range 
depending upon the fraction of non-DCC pions present in the sample. For 
various values of /3, and total number of non-DCC pions, (1 — [3)N n , the above 
equation can be evaluated easily to obtain the final / distribution from which 
one can get the resultant fluctuation, < 8f 2 >. Knowing < 5f 2 > the $- 
measure can be easily estimated. The results obtained for an average total 
number of 300 pions in an event, are shown in Fig. 1. One can see that the 
value of ^-measure, decreases with decrease in the fraction of DCC pions. 

Multiple domains - It is possible that in an event multiple domains of DCC 
are formed. In such a case the total probability distribution of the neutral 
pion fraction is the average value of P(f) over all of the domains. This can be 
written as 



P m (f) = /#!••• df m S(f ~ m + )Pl(fl) • • • Pm(fm) (15) 

where, m is the number of domains. It can be shown that the resultant prob- 
ability distribution approaches a gaussian centered at 1/3 with the standard 
deviation ~ l/y/m. This means < 5f 2 > ~ 1/m, so that $, reduces as the 
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number of domains increases in an event. 



^muiti-Dcc ~ \J< K > /m - (16) 

However, by carrying out this analysis by dividing the n — phase space to 
appropriate bins this effect can be reduced. 

Neutral pion decay. - The neutral pions decay to photons (n° — > 27), before 
they reach the detectors. For such a case one can take N n o = 2N 1 . Then 
carrying out the analysis similar that done for arriving at Eqn. 10, one obtains 

% = y/A < N n > < 5P > - v /2</>(1-</>) (17) 

This shows decay introduces a finite value of $ even for non-DCC, of the order 
of 0.27 

Detector effects - We know that the measurement of any observable in an 
experiment is affected by the efficiency of the detector. The effect of this can 
be seen easily through the following calculation. Consider the efficiency of 
detecting photons is e\ and that for charged pions is e 2 . Then we have, N* xp = 
€\N 1 and N^ p = €2N c h. Photon multiplicity measurements are also affected 
by charged particle contamination. But true iV 7 can be obtained from the 
measured iV 7 _ likc using the relation [27] 

iV 7 = J- iVlike (18) 

where, e 7 is the photon counting efficiency, p 1 is the purity of the photon 
sample obtained from detector simulations. For convenience one can define 
I 2 as ei- Then following the simple statistical analysis and assuming that the 
fluctuation in efficiencies are independent of multiplicity (for simplicity and 
may indeed be true if the analysis is carried out in narrow bins in centrality) 
and fluctuations in efficiencies are small, we find 

$e* P = 2e ie2V /< N n >< Jp >+ <N n xf > 2 < (1 -f)> 2 E (19) 

-^2e l e 2 < f > (1- < / >) 

where, E=[^ + |$]. 

Let us consider an experiment where the average total detected pion multi- 
plicity is 300. Taking typical values of relative fluctuation (a /mean) in effi- 
ciencies, to be ~ 3% [28], we have < 5e x > 2 /< e\ > 2 and < 5e 2 > 2 /< £2 > 2 
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both to be about 0.0009. Together they introduce an error 0.0027. Then the 
term < N n >< f > 2 < (1 - /) > 2 [||^ + ||^] ~ 0.04 which is very small 
compared to < N n > < 5f 2 > so may be for simplicity neglected. Hence Q exp 
can be now written as, 

§ exp = 2e lC2V /< N n >< 5P > - ^2t x e 2 < f > (1- < / >) (20) 

So knowing the efficiency and purity of the detected multiplicity sample, one 
can estimate the value of $ exp . 

Application to simulated data - We have discussed above many factors that 
possibly affect the detection of DCC-type fluctuations in data. Now lets apply 
the ^-measure to simulated data and observe the sensitivity of the measure 
to detection of DCC-type fluctuations. Another advantage of applying it to 
simulated data, is that it can form a guidance for application of this to actual 
data in future. Simulated events were generated using VENUS 4.12 [29] event 
generator with default parameters. 15K VENUS events with impact parameter 
less than 3 fm was generated for this study. We have compared the observed 
effect to simulation results based on a simple DCC model [12,30]. We assume 
the formation of a single DCC domain of a given size lying within the limited 
coverage (77 = 3-4, and full azimuthal) of a hypothetical detector system, 
consisting of a photon multiplicity detector and a charged particle multiplicity 
detector. In this model the output of the VENUS 4.12 event generator has been 
suitably modified to accommodate such a domain, characterized by percentage 
of pions being DCC type. To introduce DCC in a certain fraction of pions out 
of the total pions in an event, the charge of the pions is interchanged pairwise 
(7r + 7r~ <-> 7r°7T ), according to the DCC probability as given by Eqn. (1). After 
allowing the 7r°s to decay, all the particles are passed through the realistic 
detector responses. The efficiency of charged particle detection was taken to 
be 90% ±5% [31]. The photon counting efficiency and purity of photon sample 
was taken to be same and 70% ± 5% [27]. The average number of photons 
detected are 372 and the average number of charged particles detected are 
456. 

First we study the effect of number of pions being DCC-type. As discussed 
above and shown in Fig. 1, the value of $ decreases as the number of non- 
DCC pions in a given event increases. Here we study the same in a realistic 
scenario using simulated data. We introduce DCC for different fractions of 
pions being DCC-type ((3) in an event. This is done following the method 
discussed above. Then we calculate the quantity, $ — $ norma ;, Where <3> is the 
measure of fluctuation for a given DCC pion fraction ((3) and ^normal is the 
measure of fluctuation for normal or non-DCC type events (detector and decay 
effects included). The variation of $ — ^normal as a function of (5 is shown in 
Fig. 2. From the Fig. 2 we find that the fluctuations decreases with decrease 
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Fig. 2. Variation of — ^normal as a function of the pion fraction {0) from a DCC 
origin, obtained from simulated data using VENUS. 

with DCC pion fraction, which is consistent with the theoretical calculation 
shown in Fig. 1. The statistical error on <£> measure was calculated by taking 
different number of non-DCC type events (1000, 2000, • • -,10000) then finding 
the maximum variation in $ value for these sets. It was calculated to be 0.006. 
From this we find that such a typical experiment as discussed in the generating 
simulated data, will be able to detect DCC-type fluctuation, if the number of 
pions in an event being DCC-type is above ~ 3%. 

Secondly, we know that all events in a given heavy-ion reaction cannot be of 
DCC type. The effect of this has been studied for a fixed number of pions 
being DCC-type. For simplicity, we have assumed that for DCC type events 
the percentage of DCC-type pions out of the total number of pions is 25%. 
Then we have varied the number of events being DCC-type (a) out of the 
total number of events in an given ensemble of events, to study the effect 
of DCC event fraction (a) on ^-measure. The results have been presented as 
$ — ^normal vs - a and are shown in Fig. 3. We find that the value of $-$„ om(1 i 
decreases with decrease in a, which is along expected lines. Keeping in mind 
the above mentioned statistical error on ^-measure, we find that for DCC 
pion fraction {(3) of 0.25, the measure is sensitive to more than 1% of events 
being DCC-type out of a given ensemble of events. Similar calculations can 
be done for different values of f3, but not discussed here, as the aim here is to 
demonstrate the utility of this widely used observable for DCC-type of studies. 
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DCC event fraction (a) 



Fig. 3. Variation of $ — Q nor mai as a function of the DCC event fraction (a) for a 
fixed DCC pion fraction ((3 = 0.25), obtained from simulated data using VENUS. 



Summary - We have discussed the utility and sensitivity of ^-measure for look- 
ing at DCC-type fluctuations. We have studied the properties of ^-measure for 
three different simple models of multiplicity distributions. Then we discussed 
the effect of various factors affecting the DCC-type fluctuations on ^-measure. 
These include the existence of background or non-DC C type pions in addition 
to those from DCC origin in an event, possible existence of multiple domains 
of DCC, decay of neutral pion to photons and the detector effects like effi- 
ciency and purity. We have then discussed extensively the variation of $ with 
fraction of pions being DCC type (f3) and fraction of events being DCC type 
(a), within the framework of a simple DCC model. From the study using simu- 
lated data from VENUS and incorporating realistic detector effects, we found 
that for a typical example, such a measure is sensitive to DCC-type fluctua- 
tions with (3 greater than 0.03 and for a f3 — 0.25, it is sensitive to DCC-type 
fluctuation with a greater than 0.01. With this study we have demonstrated 
that ^-measure is a sensitive observable to look for DCC-type fluctuations, 
while it preserves its other properties, which has made it a powerful measure 
in event-by-event studies. 
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